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Abstract 

We perform a global fit to the CKM unitarity triangle using the latest experimental and theo- 
retical constraints. Our emphasis is on the hadronic weak matrix elements that enter the analysis, 
which must be computed using lattice QCD or other nonperturbative methods. Realistic lattice 
QCD calculations which include the effects of the dynamical up, down, and strange quarks are now 
available for all of the standard inputs to the global fit. We therefore present lattice averages for all 
of the necessary hadronic weak matrix elements. We attempt to account for correlations between 
lattice QCD results in a reasonable but conservative manner: whenever there are reasons to believe 
that an error is correlated between two lattice calculations, we take the degree of correlation to 
be 100%. These averages are suitable for use as inputs both in the global CKM unitarity triangle 
fit and other phenomenological analyses. In order to illustrate the impact of the lattice averages, 
we make Standard Model predictions for the parameters Bk, \V c b\, and |T4b|/|Kb|- We find a 
(2-3)<r tension in the unitarity triangle, depending upon whether we use the inclusive or exclusive 
determination of |V^,|. If we interpret the tension as a sign of new physics in either neutral kaon or 
B-mixing, we find that the scenario with new physics in kaon-mixing is preferred by present data. 
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I. INTRODUCTION 



The B-factories and the Tevatron have been remarkably successful in producing a wealth 
of data needed to constrain the flavor sector of the Standard Model. Inconsistencies between 
independent determinations of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [1, 2] and its CP- violating phase would provide evidence for new physics. The search 
for such inconsistencies has received a great deal of attention because generic new physics 
scenarios lead to additional CP-violating phases beyond the single one of the Standard 
Model. Physicists generally believe that the Standard Model cannot be the whole story, 
despite its great experimental successes, because it does not describe the large matter- 
antimatter asymmetry of the universe, nor can it account for dark matter. Although there 
is reasonably good agreement with the Standard Model prediction of a single CP-violating 
phase, as encoded in the CKM matrix, some tensions have been pointed out recently [3-7]. 

Most of the constraints on the CKM matrix are limited by theoretical uncertainties in 
the hadronic matrix elements that encode the nonperturbative QCD contributions to weak 
processes. These hadronic matrix elements are calculated using numerical lattice QCD, 
and great progress has been made in reducing the errors in the last few years. 1 Three 
flavors of light quarks (u, d, and s) are now being included in the vacuum polarization, 
so that the quenched approximation has become a thing of the past for most quantities of 
interest. Lattice calculations of many quantities have now been done by various groups with 
all sources of systematic error under control. 2 In order to maximize the impact of lattice 
input on phenomenology, it is necessary to average the different results. This is not entirely 
straightforward, since correlations between various lattice errors must be taken into account. 
For example, statistical errors in quantities computed on the same gauge ensembles will be 
highly correlated. Systematic errors can also be correlated, so familiarity with the lattice 
methods used in each calculation is needed in order to understand and account for these 
correlations in the averaging procedure. 

In this work we present lattice QCD averages of the hadronic weak matrix elements 
that enter the global fit of the CKM unitarity triangle. We provide results for the neutral 
kaon mixing parameter (Bk), for the neutral 5-meson decay constants and mixing matrix 

1 For a pedagogical review of lattice QCD methods see Ref. [8]. 

2 See, for example, Refs. [9, 10], for the status of lattice calculations of kaon and heavy-light physics. 
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elements, for the inclusive determinations of the CKM matrix elements |T4&| and \ V c b\, for the 
Standard Model correction to e K (/t e ), and for the kaon decay constant (fx)- Although we do 
not know the exact correlations between different lattice calculations of the same quantity, 
we still attempt to account for correlations in a reasonable manner. We do so by assuming 
that, whenever a source of error is at all correlated between two lattice calculations, the 
degree-of-correlation is 100%. This assumption is conservative and will lead to somewhat 
of an overestimate in the total error of the lattice averages; nevertheless, it is the most 
systematic treatment possible without knowledge of the correlation matrices between the 
various calculations, which do not exist. These averages are intended for use in the global 
CKM unitarity triangle fit, as well as other phenomenological analyses. 3 

Our emphasis in this work is different than that of the CKMfitter and UTfit collab- 
orations [11, 12]. Their focus is primarily on the statistical techniques used to extract 
information from a given set of inputs (see Ref. [13] for a quantitative comparison of the 
Bayesian versus frequentist approaches), whereas our focus is on the lattice QCD inputs 
themselves. Nevertheless, it is useful to point out some important differences between the 
manner in which they compute their lattice averages and the treatment that we use in this 
paper. Both CKMfitter and UTfit combine two- and three-flavor results in their lattice 
averages [14, 15], which we do not. This is because, although the systematic uncertainty 
in a particular quantity due to omitting the dynamical strange quark may in fact be small, 
this is impossible to quantify until the equivalent calculation has been done with both two 
and three flavors. We therefore only consider Nf = 2 + 1 lattice results in our averages. 
CKMfitter assigns the smallest systematic error of any of the individual lattice calculations 
to the average [14] instead of combining the systematic uncertainties between different lat- 
tice calculations. This treatment does not take full advantage of current lattice QCD results 
by preventing the average value from having a smaller systematic uncertainty than any of 
the individual values. Thus, although this treatment is conservative, it may obscure the 
presence of new physics in an attempt to be overly cautious. The method for obtaining the 

3 There are some additional correlations between different lattice quantities that enter the unitarity trian- 
gle fit, since collaborations often calculate more than one quantity using the same gauge configurations. 
Although we do not include these effects in our analysis, such correlations are reduced in the procedure 
of averaging several lattice results for the same quantity that have been computed with different con- 
figurations. Therefore these correlations should become less significant over time as more independent 
Nf = 2 + 1 flavor lattice results become available to include in the averages. 
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central values and errors used by UTfit is not fully spelled out in Ref. [15]. In this work 
we take all quoted lattice errors at face value when averaging results, but we only include 
results with complete systematic error budgets that have been sufficiently documented in 
either a publication or conference proceeding. 

This paper is organized as follows. In Sec. II we average lattice QCD results for hadronic 
weak matrix elements that enter the standard global fit of the CKM unitarity triangle. 
We present the individual results used in the averages, briefly describe the methods used 
in each lattice calculation, and spell out which errors we consider correlated between the 
calculations. Next, for completeness, we summarize the other inputs used in our unitarity 
triangle fit in Sec. III. We then illustrate the impact of the new lattice averages on the 
unitarity triangle fit in Sec. IV. We observe a (2-3) a tension in the fit, depending upon 
whether we use the inclusive or exclusive determination of \V c b\. In Sec. V we interpret the 
tension as a sign of new physics in either neutral kaon or B-meson mixing. We find that the 
current data prefer the scenario in which the new physics is in the kaon sector. Finally, we 
summarize our results and conclude in Sec. VI. Summary plots of all of the lattice averages 
are provided in Appendix A. 

II. LATTICE QCD INPUTS TO THE FIT OF THE UNITARITY TRIANGLE 

Many of the constraints on the CKM unitarity triangle rely upon knowledge of hadronic 
matrix elements that parameterize the nonperturbative QCD contributions to weak decays 
and mixing. In the past, these hadronic weak matrix elements have often been difficult to 
compute precisely, and have enabled only mild (10 — 15% level) constraints on the apex of 
the CKM unitarity triangle that have been insufficient to probe the presence of new physics. 
Recent advances in computers, algorithms, and actions, however, now allow reliable lattice 
QCD calculations of hadronic weak matrix elements with all sources of systematic error 
under control. State-of-the-art lattice computations now regularly include the effects of the 
dynamical up, down, and strange quarks. They also typically simulate at pion masses below 
300 MeV, and sometimes even below 200 MeV, in order to control the extrapolation to the 
physical pion mass. For many hadronic weak matrix elements of interest, there are now at 
least two reliable lattice calculations. Just as with experimental measurements, some of the 
errors are correlated among the lattice QCD calculations, and such correlations much be 
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B K 


(5 Bk) stat 


{SB^syst 


HPQCD/UKQCD '06 [17] 


0.83 


0.02 


0.18 


RBC/UKQCD '07 [18] 


0.720 


0.013 


0.037 


Aubin, Laiho & Van de Water '09 [19] 


0.724 


0.008 


0.028 


Average 


0.725 ±0.026 







TABLE I: Unquenched lattice QCD determinations of the neutral kaon mixing parameter Bk- A 
plot showing the three Nf = 2 + 1 results and their average is given in Fig. 5. 

taken into account when averaging lattice inputs to be used in the CKM unitarity triangle 
analysis. 

In this section, we average the latest lattice QCD results and provide values that should be 
used in current fits of the CKM unitarity triangle. In the averages, we only include results 
from simulations with three dynamical quark flavors, and with associated proceedings or 
publications that include comprehensive error budgets. Fortunately, for all quantities of 
interest, there is at least one calculation that satisfies these critera. In taking the averages 
we assume that all errors are normally distributed and follow the prescription outlined in 
Ref. [16] to take the correlations into account. The degree of correlation induced by a given 
source of uncertainty onto the errors of different lattice calculations is extremely difficult 
to estimate. In order to be conservative, whenever there are arguments that suggest some 
correlation between errors in distinct lattice results, we take it to be 100%. Finally, we adopt 
the PDG prescription to combine several measurements whose spread is wider than what 
expected from the quoted errors: the error on the average is increased by the square root of 
the minimum of the chi-square per degree of freedom (constructed following Ref. [16]). 

A. B K 

The experimental measurement of indirect CP- violation in the kaon sector, ek, when 
combined with a nonperturbative determination of the neutral kaon mixing parameter, Bk, 
places a constraint on the apex of the unitarity triangle. There have been three realistic 
lattice QCD calculations of B K since 2006; the results are summarized in Table I. 

The first, by the HPQCD and UKQCD Collaborations [17], uses the "2+1" flavor asqtad- 
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improved staggered gauge configurations [20-22] generated by the MILC Collaboration [8, 
23], which include the effects of two degenerate light quarks and one heavier quark with a 
mass close to that of the physical strange quark. The calculation also uses staggered valence 
quarks in the four-fermion operator used to compute Bk- The result for the renormalization- 
group invariant quantity Bk has a ~ 22% total uncertainty, which is primarily due to the 
omission of operators specific to staggered fermions that break flavor symmetry in the lattice- 
to-continuum operator matching calculation. Because the other determinations of B K have 
much smaller total errors, this result has little impact on the weighted average. 

The second calculation by the RBC and UKQCD Collaborations [18] uses 2+1 flavor 
domain-wall gauge configurations, as well as domain-wall valence quarks in the four-fermion 
operator used to compute Bk- Because domain- wall quarks have an approximate chiral 
symmetry [24, 25], it is easier to calculate the renormalization factor needed to determine 
B K in the continuum and in the MS scheme for domain-wall quarks than for staggered 
quarks. Therefore the RBC and UKQCD Collaborations compute the renormalization factor 
in the RI-MOM scheme nonperturbatively using the method of Rome-Southampton [26], and 
convert it to the MS scheme using 1-loop continuum perturbation theory [27]. They estimate 
that the truncation error due to the perturbative matching is small, ~ 2%. Although the 
RBC/UKQCD result is obtained from data at only a single lattice spacing of a m 0.11 fm, 
they include a reasonable ~ 4% estimate of discretization errors in their total uncertainty, 
which is based on the scaling behavior of quenched data with the same gluon and valence 
quark action. The total error in the RBC/UKQCD calculation of B K is ~ 5%. 

Recently, Aubin, Laiho, and Van de Water (ALV) calculated B K using domain-wall va- 
lence quarks on the MILC staggered gauge configurations [19]. The use of domain-wall 
valence quarks allows them to compute the renormalization factor nonperturbatively in 
the RI-MOM scheme just like RBC/UKQCD. Their result, however, includes a more con- 
servative estimate of the truncation error, ~ 3%, which is based on a comparison with 
an independent calculation of the renormalization factor using lattice perturbation theory. 
Their calculation also improves upon the work of RBC and UKQCD by analyzing data at 
two lattice spacings, so that they can extrapolate Bk to the continuum limit. They obtain 
a total error in B K of ~ 4%. 

In order to average the three determinations of B K , we must determine which sources 
of error are correlated between the various calculations. The HPQCD/UKQCD and ALV 
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calculations both use the staggered gauge configurations generated by the MILC Collabora- 
tion. Aubin, Laiho, and Van de Water, however, use nine independent ensembles of gauge 
configurations for their analysis, whereas HPQCD/UKQCD only use two ensembles. Thus 
the overlap in the two sets of data is quite small, and the statistical errors of the ALV result 
are sufficiently independent of the statistical errors of the HPQCD/UKQCD result that we 
treat them as uncorrelated in the average. The RBC/UKQCD and ALV calculations both 
use the same 1-loop continuum perturbation theory expression to convert the B K renormal- 
ization factor from the RI-MOM scheme to the MS scheme. Thus we treat the truncation 
errors as 100% correlated between the two calculations. Given these assumptions, we obtain 

B K = 0.725 ± 0.026 (1) 

for the weighted average, and we use this value for the unitarity triangle fit presented in 
Sec. IV. 

B. S-meson decay constants and mixing matrix elements 

The 5-meson decay constant f B places a constraint on the CKM unitarity triangle when 
combined with the experimental branching fraction for B — > tv leptonic decay. Because 
the experimental measurement is difficult, the B — > tu unitarity triangle constraint is cur- 
rently quite weak, and is not included in the standard global unitarity triangle fits [11, 12]. 
Nevertheless, we present the 5-meson decay constant here with the expectation that the 
experimental branching fraction will improve and the constraint will be more useful in the 
future. Furthermore, we can use the average values for the Bd and 5 s -meson decays con- 
stant, which have smaller errors than any of the individual determinations, to reduce the 
total uncertainty in the Bd and £> s -meson mixing matrix elements, as discussed later in this 
section. There have been two 2+1 flavor lattice calculations of the _B-meson decay constants; 
the results for f B are summarized in the upper panel of Table II, while those for f Bs in the 
lower panel. 

The Fermilab Lattice and MILC Collaborations presented preliminary determinations of 
f B and f Bs at Lattice 2008 [28]. They use the asqtad action for the light u, d, and s-quarks, 
and the Fermilab action [30] for the heavy 6-quarks. The largest uncertainty in the decay 
constants is statistical, and is ~ 3% for both f B and f Bs . Several systematic uncertainties 



7 



/ s (MeV) (5/ B )stat (Sfs) syst 

FNAL/MILC '08 [28] 195 7 9 

HPQCD '09 [29] 190 7 11 
Average 192.8 ± 9.9 

J syst 

FNAL/MILC '08 [28] 243 6 9 

HPQCD '09 [29] 231 5 14 
Average 238.8 ± 9.5 

TABLE II: Unquenched lattice QCD determinations of the £?-meson decay constants fs and fs s ■ 
Plots showing the Nf = 2 + 1 results and their averages are given in Figs. 6 and 7. 



- the light-quark discretization error and chiral extrapolation, heavy-quark discretization 
error, and scale and light-quark mass determination — all lead to comparable errors of ~ 
2%. 

The HPQCD Collaboration recently published a determination of /# and fs s [29] using 
staggered light quarks and NRQCD 6-quarks [31]. The statistical plus chiral extrapolation 
errors are comparable to those of Fermilab/MILC. The largest systematic errors, however, 
are from the continuum extrapolation (~ 3%) and operator matching (~ 4%). 

Because both decay constant calculations rely upon the MILC gauge configurations, in- 
cluding many overlapping ensembles, we treat the statistical errors as 100% correlated be- 
tween the two calculations. Most of the systematic errors in the two calculations, however, 
such as those from tuning the quark masses, heavy-quark discretization effects, and operator 
matching, are independent, so we treat the systematic errors as uncorrelated. Given these 
assumptions, we obtain the weighted averages 

f B = 192.8 ±9.9 (2) 
f Bs = 238.8 ±9.5. (3) 

In practice, the CKMfitter and UTfit Collaborations do not in fact, use the 5-meson decay 
constant to implement the unitarity triangle constraint from B — > tu decay. Instead, they 
construct the ratio B.R.(£> — > Tv)/Amd, where Am^ is the i^-meson oscillation frequency, 
to reduce the uncertainty from hadronic matrix elements. The quantity f B cancels in this 
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HPQCD '09 [29] 1.26 ±0.11 


1.33 ±0.06 



TABLE III: Unquenched lattice QCD determinations of the neutral S-meson bag parameters Bb q - 



ratio, such that the ratio depends only on the _B-meson bag parameter, Bs d , which currently 
has a smaller relative uncertainty than f B . Currently there is only one available 2+1 flavor 
calculation of the neutral 5-meson bag parameters by the HPQCD Collaboration [29]. They 
use the same lattice actions and analysis methods as for the decay constants, and obtain 

B Bd = 1.26 ±0.11 (4) 
B Bs = 1.33 ±0.06. (5) 

These results are also presented in Table III. 

The experimental measurements of the B^- and 5 s -meson oscillation frequencies, when 
combined with a calculation of the neutral 5-meson mixing matrix elements, place additional 
constraints on the apex of the CKM unitarity triangle. The weaker of the two constraints 
comes from Am^, which is proportional to the hadronic matrix element fs d \j Bs d - Nev- 
ertheless, this constraint plays an important role in the search for new physics because, 
depending upon the type of physics beyond the Standard Model (BSM) that is present, new 
physics may affect B s - and 5^-mixing independently. For example, in some minimal flavor 
violating scenarios, new physics will alter the separate constraints on the apex of the CKM 
unitarity triangle from B s - and .B^-mixing, but not the constraint from their ratio. 

Although there has been only one 2±1 flavor calculation of the neutral 5-meson mixing 
matrix elements by the HPQCD Collaboration [29], there have been two calculations of the 
decay constant as discussed earlier in this section. We can therefore use the average 
values of J'b and fs a to improve the lattice determinations of the mixing matrix elements 
fs d \J Bb a and fs s \J Bb s - We do so by combining the average value of fs in Table II 
with the HPQCD determination of Bb a in Table III. This procedure reduces the errors 
in the mixing matrix element to below that from the HPQCD calculation alone, thereby 
improving the resulting constraint on the unitarity triangle. We add the errors of fs and 
B Bd in quadrature, despite the fact that the average fs value contains information from the 
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f B ^B Bd (Me\) f Bs ^B Bs (MeV) 
Average 216 ± 15 275 ± 13 

TABLE IV: Unquenched lattice QCD averages of the neutral £>-meson mixing matrix elements 



f B y B Bd and f Ba \jB Ba . The results are obtained by combining the average decay constants given 
in Table II with the HPQCD determinations of the bag-parameters presented in Table III, thereby 
minimizing the total uncertainties. 

HPQCD decay constant calculation, and is therefore somewhat correlated with the HPQCD 
B Bd value. This error treatment is conservative, however, because adding the HPQCD 



errors for f B and B Bd in quadrature slightly overestimates the resulting error in f Bd yB Bd 
since correlated statistical fluctuations lead to a smaller error in the product of the two 
quantities (this is also true for the £?<,-meson case). The resulting 2+1 flavor lattice averages 



for f Bd \jB r B ~ d and f Bs \Jb\ s are given in Table IV. 



In practice, the hadronic matrix element fB d \jB Bd has larger uncertainties than the 



corresponding quantity in 5 s -mixing, f Bs yB Bg . This is primarily because current lattice 
QCD calculations can simulate directly at the physical s-quark mass, but must extrapolate 
to the u- and rf-quark masses. Therefore the chiral extrapolation error, which is often the 
dominant systematic, is larger for f Bd \J B Bd than for f Bs \J B Bs . In order to minimize the 
hadronic uncertainty in the Am d constraint on the unitarity triangle, we therefore replace 



f'B d \JB Bd with f Bs \JB B Ji, where f = f Bs y B Bs / f Bd y B Bd is an S'f/(3)-breaking ratio 
that can currently be determined more accurately in lattice calculations than the individual 
matrix elements, as discussed below. 

The more stringent neutral 5-meson mixing constraint on the unitarity triangle comes 
from the ratio of the oscillation frequencies, Am s / Arrid, because many uncertainties are 
reduced in the lattice calculation of £. There have been two recent 2+1 flavor lattice QCD 
calculations of £; the results are summarized in Table V. 

The Fermilab Lattice and MILC Collaborations presented a preliminary calculation of £ 
at Lattice 2008 [32], while the HPQCD Collaboration has recently published a determination 
of £ to 2.6% [29] accuracy. Both groups use the same lattice actions and methods as they 
did for the 5-meson decay constant calculations. The largest uncertainties in Fermilab- 
MILC's determination of £ are from statistics and from the chiral-continuum extrapolation, 
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£ OOstat (^)syst 

FNAL/MILC '08 [32] 1.205 0.036 0.037 

HPQCD '09 [29] 1.258 0.025 0.021 

Average 1.243 ± 0.028 

TABLE V: Unquenched lattice QCD determinations of the 5L r (3)-breaking ratio £. A plot showing 
the two Nf = 2 + 1 results and their average is given in Fig. 8. 



both of which are ~ 3%. They obtain a total error in £ of ~ 4%. HPQCD's largest 
source of uncertainty is also statistics and the chiral- continuum extrapolation, which together 
contribute ~ 2% to the total uncertainty. Because both calculations of £ rely on the MILC 
gauge configurations, in taking the average of the two results, we treat the statistical errors 
as 100% correlated. We treat the systematic errors as uncorrelated because most of the 
systematic errors are independent since they use different heavy-quark actions and operator 
renormalization methods. Given these assumptions, we obtain the weighted average 

f = 1.243 ±0.028 (6) 

for use in the unitarity triangle analysis. 



C. \V ub \ 

The CKM matrix element \V u b\ also places a constraint on the apex of the CKM unitarity 
triangle. It can be determined by combining experimental measurements of the branching 
fraction for semileptonic B — > ^Iv decay with lattice QCD calculations of the B — > -kIv form 
factor. There have been two exclusive determinations of \V U \\ based on 2+1 flavor lattice 
calculations; the results are summarized in the upper panel of Table VI. 

In 2006 the HPQCD Collaboration published the first unquenched computation of the 
B — > ixtv semileptonic form factor using asqtad valence light quarks and NRQCD valence fa- 
quarks [33] and MILC 2+1 flavor dynamical gauge configurations. The B — > tt£u form factor 
is more difficult to compute numerically than other lattice quantities such as Bk or £, and 
consequently has a larger total error. Because of the poor statistics associated with lattice 
data at nonzero momentum, the largest source of uncertainty in the HPQCD form factor 
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38.3 


0.5 
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38.6 ±1.2 







TABLE VI: Exclusive determinations of the CKM matrix elements \V u b\ and \V c b\ from unquenched 
lattice QCD calculations. Plots showing the Nf = 2 + 1 results and their averages are given in 
Figs. 9 and 10. 

calculation is the 10% statistical plus chiral extrapolation error. When the HPQCD result 
for the form factor is combined with the latest Heavy Flavor Averaging Group (HFAG) 
average for the B — > ixtv branching fraction [34], one obtains \V u b\ with a total error of 
~ 16%, only ~ 6% of which comes from the experimental uncertainty in the branching 
fraction. 

The Fermilab Lattice and MILC Collaborations recently published an improved deter- 
mination of the B — > iriu semileptonic form factor and \V u b\ using staggered light quarks 
and Fermilab 6-quarks [35]. As in the case of the HPQCD calculation, the largest source of 
uncertainty is statistics plus chiral-continuum extrapolation, which leads to a ~ 6% error in 
the form factor. Fermilab/MILC, however, extract \V u b\ in a different manner than HFAG. 
They perform a simultaneous fit to the lattice data and the 12-bin BABAR experimental 
data [36] using a fit function based on analyticity and crossing symmetry [39-42], leaving 
the relative normalization between lattice and experiment as a free parameter to be deter- 
mined in the fit. With this method they reduce the total uncertainty in \V u b\ by combining 
the lattice and experimental information in an optimal, model-independent manner. They 
obtain a total error in \V ub \ of ~ 11%. 

Because the HPQCD and Fermilab/MILC calculations both use the MILC gauge configu- 
rations, we treat their statistical errors as 100% correlated when taking the average. We also 
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treat the experimental errors in the two determinations of \V u b\ as 100% correlated. This is 
a conservative assumption because the HPQCD extraction comes from the HFAG average, 
which is obtained from many experimental measurements of the branching fraction including 
the 12-bin BABAR analysis. We treat the systematic errors in the two calculations as un- 
correlated, since that they use different actions for the heavy quarks and different methods 
for the lattice-to-continuum operator matching. Given these assumptions, we obtain 

\V ub \ = (3.42 ±0.37) x HT 3 (7) 

for the weighted average. Note that in our averaging procedure we symmetrize the HPQCD 
systematic error. For this reason the central value of the average (7) is slightly larger than 
both the HPQCD and Fermilab/MILC central values. 

D. \V cb \ 

The CKM matrix element \V cb \ normalizes the base of the CKM unitarity triangle. There- 
fore it implicitly enters many of the constraints on the apex of the CKM unitarity triangle, 
including those coming from Bk and \V u b\. \V c b\ can be determined by combining experi- 
mental measurements of the branching fractions for B — > Div or B — > D*iv semileptonic 
decay, in combination with lattice QCD calculations of the relevant form factor at zero 
recoil. There have been two exclusive determinations of |V^,| based on 2+1 flavor lattice 
calculations; the results are summarized in the lower panel of Table VI. 

The Fermilab Lattice and MILC Collaborations presented the first unquenched lattice 
determination of the B — > Div form factor at Lattice 2004 [37]. Like other Fermilab/MILC 
calculations of heavy-light meson quantities, this calculation uses staggered light quarks 
and Fermilab b- and c-quarks. Because the B — > Div form factor at zero recoil can be 
obtained by a carefully constructed double ratio of matrix elements, it can be computed 
very precisely from lattice calculations. The statistical error in the resulting form factor 
determination is ~ 2%, and all systematic errors are ~ 1% or less. Although the result is 
obtained from data at only a single lattice spacing, Fermilab/MILC include an estimate of 
discretization effects in their error budget. When the result for the B — > Div form factor is 
combined with the latest HFAG average of the experimental branching fraction [34] , it leads 
to a determination of [V^l with a ~ 4% total error, of which ~ 2% is lattice theoretical and 
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~ 4% is experimental. 

More recently, the Fermilab Lattice and MILC Collaborations published the first un- 
quenched lattice determination of the B — > D*£u form factor [38]. This calculation also uses 
staggered light quarks and Fermilab b- and c-quarks, and obtains the form factor at zero 
recoil from a double ratio of matrix elements. It improves upon the earlier determination of 
the B — > Dtv form factor, however, by using data at three lattice spacings, and performing 
a more sophisticated chiral-continuum extrapolation. When the result for the B — > D*£u 
form factor is combined with the latest HFAG average of the experimental branching frac- 
tion [34], it leads to a determination of [V^l with a ~ 3% total error, of which ~ 2.5% is 
lattice theoretical and ~ 1.5% is experimental. We note that the average of B — > D* data 
yields a poor chi-square per degree of freedom {x^iJ d.o.f . = 39/21). For this reason, we 
rescale the experimental error quoted in Table VI by a factor ~ 1.4. 

Because the Fermilab /MILC calculations are computed on some of the same ensembles 
and using the same lattice actions and methods, we treat the theoretical errors as 100% 
correlated between the extractions of \V cb \ from D and D* final states. We assume that the 
experimental errors are independent between the two measurements. Given these assump- 
tions we obtain the weighted average 

|V c6 |exd = (38.6 ± 1.2) x 10- 3 (8) 

for use in the unitarity triangle fit. 

E. K £ 

Buras and Guadagnoli [5] have pointed out that corrections to Ek that had typically been 
neglected in the unitarity triangle analysis due to the large errors on Bk and \V c b\ (with the 
exception of Refs. [43, 44]) are actually substantial, and amount to a ~ 8% correction to the 
Standard Model prediction for e K . In this section, we provide the first estimate of the cor- 
rection factor k £ using a 2+1 flavor lattice QCD calculation of Im [A(K — > = 2))] [45]. 
We follow the notation of Ref. [5]. Our result supersedes previous estimates of k £ [46-49], 
and is in perfect agreement with the recent determination of Ref. [5] . 

It is conventional to define K — > im matrix elements in terms of definite isospin ampli- 



14 



tudes by 



A(K° -> 7T7r(/)) = A 7 e i<5 ', 



(9) 



A(a -> 7T7r(/)) = -A*e iSl . (10) 
CP-violation in the kaon system is then parameterized in terms of 



and 



2g*(<52-'5o) 

e^=— ^— w[ft-P ], (12) 



where 



_ ImAo _ Im^2 _ Re^2 

P ° = ReV P2= Re^' " = R^V (13) 

The first term in parenthesis in Eq. (11) is the short distance contribution to kaon mixing; 

it is the part that is conventionally normalized by the bag-parameter B K . The second 

term, P Q , is due to long distance contributions to kaon mixing, and is related to the ratio 

of K — > 7T7T decay amplitudes in the AI = 1/2 channel defined in Eq. (13). In the usual 

unitarity triangle analysis, <p £ is taken to be 7r/4, and P is taken to be negligible compared 

to the first term in parenthesis in Eq. (11). However, these corrections are not negligible in 

the current analysis, and the corrections coming from <f) £ ^ ir/A and Pq ^ are small but in 

the same direction [5]. 

Following Ref. [5], we define an overall multiplicative correction factor for e K that ac- 
counts for 4> £ 7^ 7r/4 and P ^ 0, 

K e = V2sin(f) e K E , (14) 

where K e parameterizes the correction to \sk\ coming from P . To good approximation, we 
have 

«« ~ 1 + -^-r» ( 15 ) 



and given that 



Re(e' K /e K )~—^(P 2 -P ), (16) 
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2+1 Flavor lmA 2 x 10 13 GeV 

RBC/UKQCD '08 [45] -7.9 ± 1.6 ± 3.9 

TABLE VII: 2+1 flavor lattice value for Im^. Errors are statistical and systematic, respectively. 
A plot comparing the Nf = 2 + 1 result with several quenched determinations is given in Fig. 11. 

we see that 

K e = 1 - -Re(e' K /e K ) + (17) 

All of the quantities entering n e , Eqs. (14) and (17), are well determined from experiment, 
assuming the Standard Model, except for P 2 . In Eq. (13), we need theory input for lmA 2 
to determine P 2 . The amplitude lmA 2 is non- vanishing due to the electroweak penguin 
contribution to K — > irir decays. We use the result of lattice calculations of this amplitude in 
our analysis. There is only one 2+1 flavor result for this quantity with rather large systematic 
errors associated with the use of leading order chiral perturbation theory [45] , obtained by 
the RBC/UKQCD Collaborations. Their value is given in Table VII, where the first error is 
statistical and the second is due to the systematic error in the determination of the leading 
order low energy constant of chiral perturbation theory (xPT), and the truncation of x?T to 
tree-level. Although the quoted error budget does not include typical lattice errors such as 
finite volume effects, scale setting, and discretization errors, they are almost certainly much 
smaller than the errors attributed to the use of leading order chiral perturbation theory. 
For comparison, we also mention lattice results for Iim4 2 in the quenched approximation, 
collected in Table VIII. Some of these quenched calculations are more thorough than others 
at assessing systematic errors. However, all of their systematic error budgets are necessarily 
incomplete, due to the uncontrolled nature of the quenched approximation, and we do not 
attempt to estimate this error here. We take a simple average of the five quenched lattice 
results and note that it is very similar to the 2 + 1 flavor result. We use only the 2 + 1 flavor 
result for our determination of K e for use in the unitarity triangle fits. 

All inputs used to determine n £ are given in Table IX. We take the most recent exper- 
imental world average for Re(e'/e) [55], noting that we innate the errors according to the 
PDG prescription because of the somewhat low confidence level (13%) in the world average. 
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Quenched lmA 2 x 10 13 GeV 



RBC '01 [50] -12.6 

CP-PACS '01 [51] -9.1 

SPQcdR '04 [52] -5.5 

Babich et al '06 [53] -9.2 

Yamazaki '08 [54] -11.8 

Average —9.6 



TABLE VIII: Quenched lattice values for ImA 2 . 



<pe = (43.51 ± 0.05)° 
to = 0.0450 

ReA 2 = 1.50 x 10~ 8 GeV 

TABLE IX: Inputs used to determine k £ . 

Using these values in Eqs. (13), (14), and (17) we find 

k £ = 0.92 ±0.01, (18) 

in agreement with Ref. [5] . The 50% error in the 2+1 flavor determination of lmA 2 dominates 
the error in K e . We note for comparison that if we use the average quenched value of Inn42, 
assigning to it a 100% error, we find K e = 0.92 ± 0.02. 

F. fx 

The kaon decay constant fx enters the CKM unitarity triangle through e K . Because 
experiments can only measure the product fx x \V US \, lattice calculations are needed to 
obtain fx by itself. There have already been four 2+1 flavor lattice QCD determinations of 
fx using different valence and sea quark actions, and several more calculations are underway. 
Thus fx is one of the best-known hadronic weak matrix elements. Table X summarizes the 
current status of 2+1 flavor lattice QCD calculations of fx- 

The MILC Collaboration published the first 2+1 flavor determination of fx in 2004 [60], 
and updated the result at Lattice 2007 by including data with lighter quarks and finer lattice 
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\e K \ = (2.229 + 0.012) x 10~ 3 
Re(e' K /e K ) = 1.68 ± 0.19 x 10~ 3 
lmA 2 = (-7.9 ± 4.2) x 10~ 13 GeV 



MMeV) (6f K ) stat (Sf K ) sys t 

MILC '07 [56] 156.5 0.4 t l 2 ° 7 

HPQCD/UKQCD '07 [57] 157 1 2 

RBC/UKQCD '08 [58] 149.6 3.6 6.3 

Aubin, Laiho & Van de Water '08 [59] 153.9 1.7 4.4 
Average 155.8 ±1.7 

TABLE X: Unquenched lattice QCD determinations of the kaon decay constant fx- A plot showing 
the four Nj = 2 + 1 results and their average is given in Fig. 12. 



spacings [56]. The largest source of uncertainty in their calculation is from the extrapolation 
to the physical light quark masses and the continuum. A small but non-negligible error also 
arises due to the determination of the absolute lattice scale needed to convert dimensionful 
quantities into physical units. MILC first determines the relative scale r\/a from the heavy- 
quark potential. Next they obtain the absolute scale r\ = 0.3108(15) (+79) by tuning f n to be 
equal to the experimental value. The uncertainty in r\ leads to an uncertainty in fx of 1"° 75, 
which is 25% of the total error. The remaining finite volume effects and EM effects are an 
order of magnitude or more smaller, and the total uncertainty in the MILC Collaboration's 
determination of f K is ~ 2%. 

The HPQCD Collaboration published a determination of f K using a mixed-action method 
with highly-improved staggered quarks [61] on the MILC asqtad-improved staggered gauge 
configurations [57]. The largest source of uncertainty in their calculation is from the deter- 
mination of the scale. They use the MILC Collaboration's determination of the relative scale 
ri/a to convert dimensionful quantities from lattice units to r\ units. They obtain the value 
T\ = 0.321(5), independently, however, using the T spectrum computed with nonrelativistic 
6-quarks on the MILC ensembles [62]. The uncertainty in r\ leads to an uncertainty in fx 
of 1.1%. The remaining statistical and systematic errors are all much smaller, and the total 
error in fx is 1.3%. 

The RBC and UKQCD Collaborations published an independent determination of fx 
using domain- wall quarks [58]. Because they obtain their result using only a single lattice 
spacing of a ~ 0.11 fm, the dominant uncertainty in their result is from discretization 
errors. They estimate these errors to be 6% using power-counting arguments. Because the 
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remaining statistical and systematic errors are all much smaller, the total error in fx is 
6.3%. 

Aubin, Laiho, and Van de Water presented a preliminary determination of fx using 
a mixed-action method with domain-wall valence quarks on the MILC staggered gauge 
configurations at Lattice 2008 [59]. The largest source of uncertainty in their calculation is 
from the chiral and continuum extrapolation, which they estimate to be 2.3%. The 1.6% 
error from the uncertainty in the scale r\ however, is close in size, so the total error in f K is 
3.0%. 

Because the HPQCD, MILC and ALV calculations all use the MILC gauge configura- 
tions, we treat the statistical errors as 100% correlated when taking the average. Because 
ALV use the MILC Collaboration's determination of the scale r\ from f n , the scale uncer- 
tainties are also 100% correlated between the calculations. We take the scale uncertainty in 
HPQCD's calculation to be uncorrelated, however, because they use a largely independent 
determination of r 1 based on the T spectrum. We also treat the remaining systematic er- 
rors as uncorrelated between the HPQCD, MILC, and ALV calculations because they use 
different valence quark formulations. The calculation of the RBC and UKQCD Collabora- 
tions is independent of the other results, and we therefore take the errors to be completely 
uncorrelated. Given these assumptions, we obtain the weighted average 

f K = (155.8 ± 1.7) MeV (19) 

to be used in the unitarity triangle analysis presented in Sec. IV. 

III. OTHER INPUTS TO THE FIT OF THE UNITARITY TRIANGLE 

Table XI summarizes the set of inputs that we use in the fit. We obtain a from the 
isospin analysis of B — > (nn, pp, pn) decays (the description of the method we use can be 
found in Refs. [63-65] and the experimental inputs are taken from Ref. [66]). We take 
the direct determination of 7 from the model-independent UTfit analysis of B — > D^K^ 
decays [67, 68] (the experimental inputs used are taken from Ref. [66]). The inclusive 
determination of \V c b\ deviates by more than 2 o from the average of exclusive results that 
we quote in Eq. (8). For use in the unitarity triangle fit we combine the three determinations 
of \V cb \ from inclusive and exclusive (D and D*) modes. Taking into account the correlations 
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iKblind = (41.31 ± 0.76) x 10~ 3 [34] |V u6 | incl = (40.3 ± 1.5±l;g) x 10~ 4 [34] 

Am Bd = (0.507 ± 0.005) ps" 1 [71] Am Bs = (17.77 ± 0.10 ± 0.07) ps^ 1 [72] 

a = (89.5 ±4.3)° 7 = (78 ± 12)° [67, 68] 

m = 1.51 ± 0.24 [73] m t:Pole = (172.4 ± 1.2) GeV [74] 

m = 0.5765 ± 0.0065 [75] m c (m c ) = (1.268 ± 0.009) GeV [76] 

7/3 = 0.47 ± 0.04 [77] e K = (2.229 ± 0.012) x 10~ 3 [70] 

T) B = 0.551 ± 0.007 [78] A = 0.2255 ± 0.0007 [79] 

S^ Ks = 0.672 ± 0.024 [66] 

TABLE XI: Inputs used in the unitarity triangle fit. Note that the most precise determination of 
m c is obtained from lattice QCD [76] . 



between the errors of the two exclusive determinations of \V c f,\ and assuming no correlation 
between inclusive and exclusive analyses, we obtain: 

l^| cxcl+incl = (40.3 ± 1.0) xKT 3 , (20) 

where the error has been appropriately rescaled following the PDG prescription. We quote 
the inclusive determination of \V u b\ from the most recent GGOU analysis [34, 69]. Be- 
cause, however, the extraction of |14&|i nc i depends strongly on the theoretical framework 
adopted [34], we adopt a conservative stance and omit |I4b|inci from the set of measurements 
that we include in the full unitarity triangle fit. Our predictions for the Standard Model 
parameters in the following section are independent of \V u b\, and our conclusions regarding 
indications of new physics in Sec. V are relatively insensitive to the value of \V u b\. Apart 
from the inputs listed in Table XI, we take Gf, ™>k, ™>w, m B d and ms s from the Particle 
Data Group [70]. 



IV. STANDARD MODEL PREDICTIONS 

In this section we extract the Standard Model predictions for Bk, \V c b\ and iKtfc/Kbl- 
We use only the three constraints from S^k s , AM Bs /AM B(j and e K) and do not include 
the constraints from \V u b\, a and 7 in the fit because predictions are almost completely 
insensitive to their impact. The analytical formulae for e K and AM B J AM Bd can be found, 
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FIG. 1: Impact of ek on the UT fit. The solid, dashed and dotted contours are obtained by 
omitting ek, S^k and AM^/AM^, respectively. The left and right panels use the exclusive and 
inclusive \V c b\ determinations, respectively. 



for instance, in Ref. [7]. 

We obtain the prediction for Bx by excluding the direct lattice determination of B^ from 
the chi-square. The dominant source of uncertainty in the extraction of Bx stems from the 
strong dependence of ek on \V c b\ [ek oc |Ki,| 4 ). This issue is even more problematic because 
of the discrepancy between the extraction of \V c b\ from exclusive and inclusive decays. For 
this reason we perform the analysis both with and without the inclusive determination of 
| V c b | . Note that when combining the inclusive and exclusive extractions of | V^, | we follow the 
PDG prescription for inflating the error when combining inconsistent measurements. We 
find: 



Km, 



'1.09 ±0.12 |v; 6 | excl 

0.903 ±0.086 |V c6 | incl ( 21 ) 
0.98 ±0.10 |K b | excl+incl 

The comparison of these predictions with the lattice determination of Bk given in Eq. (1) 
yields a deviation at the 2.9<r, 2a and 2. 4a level, respectively. We obtain the prediction 
for \V c b\ in a similar fashion by excluding the inclusive and exclusive determinations of \V c b\ 
from the chi-square. We find: 

\V cb \ &t = (43.0 ± 0.9) x 10~ 3 . (22) 

This prediction deviates by 3. Oct and 1.3<r from the exclusive and inclusive determinations 
of I V c b I , respectively. 
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Figure 1 illustrates the (2 - 3)<r tension in the fit to the unitarity triangle. In the left and 
right panels we use \V c b\ from exclusive and inclusive semileptonic B decays, respectively 
The solid, dashed and dotted contours are obtained by omitting e K , S^ K and AM B J AM Bd , 
respectively. They correspond to three scenarios in which new physics affects K mixing, the 
phase and the amplitude of Bd mixing. An alternative measure of the tension in the UT fit 
is the minimum x 2 P er degree of freedom: when we include all three constraints we obtain 
Xmin/dof = 6.1 (2.6) using | K&|cxci (inci), corresponding to a confidence level of 0.2% (7.4%). 

It is interesting to note that the errors on the fitted values of B K (~ 10% - 12%) are 
much larger than the corresponding lattice uncertainty (~ 3.5%); therefore, improvements 
on the latter are not poised to have a sizable effect on the UT tension. In contrast, the 
errors on the direct and indirect determinations of \V c b\ are similar (about 2%), indicating 
that improvements on the theoretical predictions for exclusive and inclusive semileptonic 
B decays will have a huge impact on our understanding of this issue. This discussion 
is summarized in Fig. 2, which shows the relative impact of the present |V c b| ex ci and B K 
uncertainties on the total Sk error band. 

Finally, we note that the SM prediction for the ratio is in good agreement with 

the lattice expectation and deviates by only 1.6a from the inclusive ratio: 



lib 



v.. 



cb 



0.0846 ± 0.0035 fit 

0.089 ±0.010 exclusive • (23) 
0.0969 ± 0.0068 inclusive 



V. INTERPRETATION AS NEW PHYSICS 

In this section we assume that physics beyond the Standard Model does not affect tree- 
level processes at the current level of precision, and that any sign of new physics must arise 
due to higher-order loop effects. Given these assumptions, it is well known [3-6] that the 
~ 2(7 tension in the fit to the unitarity triangle can be interpreted as a manifestation of 
new physics effects in e K and/or B d mixing. In order to test the consistency of these two 
hypotheses with the current measurements, we describe the two new physics possibilities 
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-1.0 -0.5 0.0 0.5 1.0 

P 

FIG. 2: Impact of \V c b\ (solid red line) and Bk (dashed green line) on the £k error band. The 
uncertainties induced by variations of m c , mt, r]i and k £ have a negligible impact on the £k error 
budget. 

using the following model-independent parametrization: 

£ k = C £ (e k ) sm , (24) 
Mf 2 = rl^(M? 2 ) SM , (25) 

where Mf 2 is the matrix element of the complete effective Hamiltonian between B° and B 
states. A value of C £ ^ 1 would move the location of the Ek band, while the presence of a 
Td 7^ 1 and a non-vanishing 0^ would alter the following three unitarity triangle constraints: 

AM Bd = AMgf r\ , (26) 
Pes = P + e d , (27) 
a cS = a-9 d , (28) 

where j3 e s and a e s are the angles extracted from the CP asymmetries in B —>■ J/ipK and 
B — > (7T7T, pp, pre), respectively. Although the presence of new physics effects in B s mixing 
is a very interesting possibility, we do not consider it in the model-independent analysis 
presented in this work for the following reason. Because of the smallness of the phase of the 
CKM element Vt s , any evidence of CP violation in B s decays translates immediately into 
evidence for physics beyond the SM. The golden mode studied at CDF and DO is the time 
dependent CP asymmetry in the decay B s — > J/ip<f), which presently deviates by about 2a 
from the SM expectation (~ 0). Unfortunately, however, a new phase in B s mixing does not 
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p 



FIG. 3: Fit to the unitarity triangle in the SM. We average the inclusive and exclusive determi- 
nations of \Vcb\, but use only the exclusive determination of \V u b\. The black contour is obtained 
from the minimization of the complete chi-squared. 



affect any of the observables that enter the unitarity triangle fit. Therefore, signs of new 
physics in B^/K and B s mixing can be connected only in specific new physics scenarios in 
which new phases in K, Bd and B s mixing have a common origin (see for instance Refs. [5, 7] 
for a discussion of this point). 

From the inspection of Fig. 1 we see that, as long as we consider only Ek, AMg^AM^, 
S^k s an d | Kb |) the data are not able distinguish between scenarios with new physics in 
the K or B^ sectors. In fact, both the solid (new physics in K mixing) and dashed (new 
physics in Bd mixing) contours in Fig. 1 have Xmm — 0- This tie, however, is broken by 
the inclusion of constraints on |K&| (from exclusive semileptonic b — > u decays), a (from 
B — > (7T7T, pp, pn) decays) and 7 (from B — > D^K^*' decays). Figure 3 shows the resulting 
full fit to the unitarity triangle using the combined inclusive + exclusive determination of 

|Kb|- 

In order to test the hypothesis that new physics only affects neutral kaon-mixing, we 
minimize the chi-square while excluding ek from the fit. The solid contour in the left panel 
of Fig. 4 shows the allowed (p, rf) region in this scenario. Adopting the parametrization in 
Eq. (24) we obtain the following value for the new physics contribution to ek- 





1.47 ±0.17 


1Kb 


excl ' 


(C £ ht = < 


1.21 ±0.11 


|K 6 


incl ' 




1.32 ±0.14 


1Kb 


cxcl+incl 



In the upper right and lower panels of Fig. 4 we consider scenarios in which only new physics 
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in B d mixing is allowed. For the sake of simplicity we consider only two extreme cases in 
which we take (9 d ^ 0,r d = 1) and (9 d = 0,r d ^ 1). In the former case S^ K and the 
extraction of a are affected by new physics contributions and must be excluded from the 
fit; in the latter one, only AMbJ AMB d receives contributions. The fitted values of the new 
physics parameters 6 d and r d are: 





f 

(-4.3 ±2.1)° 


r<*> cxcl 




(0.940 ± 0.036) 


rcb'excl 


(0d)fit = < 


(-2.8 ± 1.9)° 


l^cftlincl 


and (r d )st = < 


(0.950 ± 0.036) 


l^| incl (30) 




(-3.4 ±2.0)° 


Wd>\ excl+incl 




(0.946 ± 0.036) 


Wd>\ excl+incl ■ 



In this case, the tension between AM B J AM Bd , Ek and |T4&| reduces the quality of the fit: 
the fit omitting the constraint from e K has a confidence level of 91%, while the fit omitting 
the constraints from S^k and a and from AM Bs /AM Bd have a confidence level of 23% and 
30%, respectively. Thus the scenario with new physics in K mixing is favored by present 
data. This can also be seen from the inspection of Eqs. (29) and (30): C £ deviates from 
qSM _ y a t a higher confidence level than 9 d from 9 d M = 0. 

Finally, it should be noted that the marked preference for new physics in the kaon sector 
is a direct consequence of our inclusion of the lower determination of \V u b\ from exclusive 
semileptonic decays in the fit. As can be seen from the upper right-hand plot in Fig. 4, further 
removal of the iKb/K&l constraint results in a fit with a high confidence level (CL=81%). 
The overlap of the constraints from Ek and AM s /AM d , however, corresponds to a very 
large value of \V u b/V c b\ = 0.120 ± 0.017. The correlation between possible new physics in 
B d mixing and a large implied value of |V^&/Vd,| is well known, and has been discussed in 
Refs. [5, 80, 81]. 

VI. CONCLUSIONS 

Lattice QCD calculations that include the effects of the dynamical up, down, and strange 
quarks are becoming standard, and allow reliable calculations of hadronic weak matrix 
elements with all sources of uncertainty under control. Because there are now multiple 
lattice calculations of most of the hadronic matrix elements that enter the unitarity triangle 
fit, it is essential to average these results in order to reduce the theoretical uncertainties 
and obtain the most sensitive test of new physics in the flavor sector possible. We have 
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FIG. 4: Full fit to the unitarity triangle. Upper left panel: The black contour is obtained without 
the inclusion of the ek constraint. Upper right panel: The black contour is obtained without the 
inclusion of the a and (3 constraints. Lower panel: The black contour is obtained without the 
inclusion of the AM^ and AMb s / AMs d constraints. 



therefore presented averages for the hadronic weak matrix elements that enter the standard 
global fit of the CKM unitarity triangle. Although we do not know the precise correlations 
between different lattice calculations of the same quantity, we have accounted for correlations 
between the different lattice results in a conservative manner in our averages. Whenever 
there is any correlation between the statistical or a particular systematic error in different 
lattice calculations, we assume that the degree of correlation is 100%. Our lattice averages 
of hadronic weak matrix elements are therefore appropriate for use in phenomenological 
analyses such as the global CKM unitarity triangle fit. 

When these up-to-date lattice averages of the hadronic weak matrix elements are used in 
a global fit of the CKM unitarity triangle, we find a (2-3) a tension. As was first pointed out 
by Lunghi and Soni [4], this tension is primarily between the three most precise constraints 
on the unitarity triangle from sin(2/5), AM^/AMb^, and Ek, and is largely independent 
of the value of \V u b\, which differs significantly between determinations using inclusive and 
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exclusive semileptonic decays. We confirm their observation and put it on an even stronger 
footing by using lattice averages that include more recent lattice calculations and take into 
account correlations. The significance of the tension depends upon whether we use the 
exclusive or inclusive determination of \V c b\, which disagree by ~ 2a. If we assume that new 
physics does not affect tree-level processes at the current level of precision, this tension can 
be interpreted as a sign of new physics either in neutral kaon mixing or in neutral 5-meson 
mixing. We find that the current data prefer the scenario in which the new physics is in 
kaon mixing; this can be seen by the fact that the confidence level of the global fit increases 
significantly when we remove the constraint from the €k band leaving all others unchanged. 
The tension between the ex band and the other constraints is enhanced by our inclusion of 
the correction factor k € , which lowers the Standard Model prediction for ex by 8%. This 
factor has been recently included by the UTfit collaboration [82] (they find a similar tension 
in the fit), but not yet by the CKMfitter group [11]. 

The errors in the hadronic weak matrix elements needed as inputs to the unitarity triangle 
analysis will continue to decrease over the next few years, as the results included in these 
averages are updated and as new independent results using different lattice actions from 
other collaborations appear. If the tension observed in the current global fit persists as the 
theoretical errors are reduced, this may indeed be a sign of new physics. This will be difficult 
to ascertain conclusively, however, unless the inclusive and exclusive determinations of |V^,| 
converge. Thus a better understanding of the theoretical errors in both determinations is a 
high priority for flavor physics. Lattice QCD calculations of weak matrix elements are truly 
living up to their promise and may ultimately lead to the discovery of new physics in the 
quark flavor sector. 
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FIG. 5: Unquenched lattice average of the neutral kaon mixing parameter B^. The three Nj = 2+1 
lattice inputs are given by shaded green circles with error bars, while the resulting average is denoted 
as a purple vertical line with a blue error band. 



1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 


— • 


1 

• FNAL/MILC'08 
1 • HPQCD'09 


1 




• — 

1 1 1 1 1 1 1 1 1 





175 180 185 190 195 200 205 210 215 220 

f B (MeV) 

FIG. 6: Unquenched lattice average of the leptonic decay constant fg. The two Nf = 2 + 1 lattice 
inputs are given by shaded green circles with error bars, while the resulting average is denoted as 
a purple vertical line with a blue error band. 



APPENDIX A: SUMMARY PLOTS OF LATTICE QCD AVERAGES 

In this section we provide summary plots of the lattice QCD averages discussed in Sec. II. 
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